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Abstract 

The results of partial wave analysis of the aj7r~7r°-system, produced 
in the reaction 7v~Be — > 7r+27r~27r'^_Be at beam momentum of 37 GeV 
are presented. An indication of the resonant type structure with mass 
M 1.74 Gel/ and J-^Af = 0"0+ in the a;p"-system was found. The 
decay branching ratio of 02 ( 1320) ~ to ujtt~tt'' was measured. A wave 
with J^M'' = 2+1+ shows a broad bump at mass M ^ 1.7 GeV. The 
decays 7r2(1670)~ — » Lop~ and a4(2040)~ ujp~ were observed for the 
first time. The bump in the 61 (1235)7r wave with exotic quantum numbers 
J^^ = 1"+ at M ^ 1.7 GeV was observed. 



^D.V.Amelin, E.B.Berdnikov, G.V.Borisov, V.A.Dorofeev, R.I.Dzhelyadin, A.V.Ekimov, 
Yu.P.Gouz, I.A.Kachaev, A.N.Karyukhin, Yu.A.Khokhov, A.K.Konoplyannikov, 
S.V.Kopikov, V.V.Kostyukhin, S.A.Likhoded, V.D.Matveev, A.P.Ostankov, D.I.Ryabchikov, 
A.A.Solodkov, O.V.Solovianov, E.A.Starchenko, N.K.Vishnevsky, E.V.VIasov, A.M.Zaitsev 
Institute for High Energy Physics, 142284, Protvino, Russia 
G.Sekhniaidze, E.G.Tskhadadze 
Institute of Physics, 380077, Tbilisi, Georgia 



1 



1 Introduction 



The results of partial wave analysis (PWA) of the wTr^Tr'^-system, produced in 
the reaction ■k~ Be — > ujTT~Tr'^Be by pion beam with momentum P^r — 37 GeV 
are presented in this article. These results were partially reported at the con- 
ferences §, [|. 

The study of this reaction is an interesting subject for several reasons. The 
quark model predicts many meson resonances with masses in the region of 1.5 — 
2.5 GeV. The exotic mesons with the same masses may exist apart from the 
ordinary qq states. The hybrids, four-quark mesons, hadron molecules can be 
in the isovector state. Practically all of them should decay to a;7r~7r°, where 
the wp", 6i7r, p3(1670)7r should be the main chanels. There exist in several 
cases definite predictions for probabilities of the decay of different objects into 
these chanels, namely wp". Study of these decays can play important role in 
determination of the nature of states under study, e.g. the results of the model 
calculations for decay widths of the 7r(1800) (|],@. That state can be a hybrid 
or a radial excitation of the tt- meson, S^Sq. The difference in the decay widths 
into ujp^ ^ predicted by the ^Pq model makes it possible to distinguish between 
these alternatives. Similar predictions were made ||^ for the ai(1700) state 
which can be a hybrid or a 2^Pi radial excitation state. The predictions for 
the 7r2(1670) and for more heavy states with J^M"^ = 2^0+ ||] are worth 
mentioning. 

The system composed of two vector mesons has been previously studied 
in the 7 — 7 interactions j6|, in the J/^ radiative decays |Q, in the nucleon- 
antinucleon annihilation , in the charge exchange processes , in the central 
[ p^ and diffractive production jlj. The relatively small probability of the pro- 
duction processes and/or the high multiplicity of the final states did not let them 
to take the experimental data sample sufficient for making detailed PWA. Nev- 
ertheless, the indications to the existence of isoscalar resonances with quantum 
numbers 2++ in the pp and ujuj systems were obtained. 

The experimental information about isovector mesons is even scarcer. The 
diffractive productions of the isovector resonances with = 0~, l"*", 2~... is 
especially suitable, for the production cross sections are large and the physical 
backgrounds are small. The waves with quantum numbers — 2+, 4+ become 
important at relatively large momentum transfer. 

2 Set-up description 

The measurements were carried out using VES spectrometer exposed to the 
negatively charged pions with momentum 37 GeV. VES is a wide aperture 
spectrometer of charged particles and photons. The beam part is composed of 
the threshold Cerenkov counters and a system of propotional chambers. The 
vertex part consists of a 4cm thick Be target and a guard system surrounding it 
to veto events with particles moving backward to the beam in the event centre 
of mass system. The spectrometer of charged particles consists of a tracking 
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system, a wide aperture magnet( kick momentum 0.7 GeV ) and a multichanel 
threshold Cerenkov counter. A fine grained electromagnetic calorimeter (1600 
cells) measures 7 energy and position. A triggering system, which consists of 
a hodoscope, scincillation counters and a guard system selects events of beam 
interaction with the target and at least two tracks in the forward spectrom- 
eter. The data sample used for presented analysis included ^ 6 ■ 10^ recon- 
structed events, corresponding to the integrated luminosity of ~ 100 fib~^. The 
acceptance for events with several charged and neutral pions is rather high { 
e « 0.6 for Li;(7r+7r^7r*')7r^7r° at M^j,-^a = l.bGeV ) and is a smooth function 
of kinematical variables. The measured width of the tt" meson in the 77 in- 
variant mass spectrum charactirises the experimental resolution of calorimeter 
and is equal to ~ 20 MeV. The measured width of the w-meson F ~ 30 MeV 
includes the resolution of the tracking part of set-up. 

3 Event selection and main features of the reac- 
tion TT^Be UTT^TT^Be 

At first, the events of the process ir^Be ■k'^2'k~2'k^ Be were selected for fur- 
ther analysis of the Lon~-K^ system using the following criteria for reconstructed 
events: 

• Two negative and a positive track originated in the target region. 

• At least four energy clusters in the calorimeter not associated with tracks. 
A cluster is assumed to be a 7-shower, if it satisfied the criterion on the 
transverse form of the shower. Events with the pairs of photons with 
invariant mass nearest to the PDG value and within 60 MeV range for tt" 
and 156 MeV for r] were selected. 

• Some ~ 10% of events with electrons were rejected using the ratio of the 
energy deposited in calorimeter and the momentum measured by magnetic 
spectrometer. 

• The Cerenkov counter responce to passed track was applied for charged 
K-mesoTL identification. Some 3.4% of events were rejected as having 
kaons. 

• Events with two tt" and, possibly free soft photons {E^ < 1 GeV), and the 
total forward final state particles energy in the range 34 GeV <^Ei < 
?>%GeV were classified as those of type tt~ Be 7r+27r~27r°i3e . 

A IC-fit was done on the mass of tt^ . The invariant mass distribution of the 
Tr+TT^vr*' system with additionaly produced tt^tt" is presented in fig. |l|. Clear 
peaks corresponding to the 77 and w-meson production with large backgound 
are evident. The a>-meson peak has a non-gaussian form. Its parameters were 
determined by the fit with the sum of two gaussian functions in the range from 
0.6 GeV to 0.9 GeV. The mass of the narrow part was M^ ^ 783.37±0.06 MeV, 
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where the systematic error was estimated to be gm^ — 1 MeV. The width was: 
Ti^ = 26.2 ± 0.2 MeV and the number of events with uj production turned out 
to be (361 ± 1) • 10^. The fraction of background events at the point of the peak 
maximum was estimated to be 30%. 

To make various distributions describing the production and decay of the 
ijj'K~'K^ system the background was estimated by the neighbouring to the uj 
peak bands and was subtracted. Apphcation of this method assumes smooth 
behaiviour of the background. The drawback is the difference of the phase 
spaces for those bands, where the background was estimated, and for the signal 
ones. 

Another method of background subtraction was apphed for cross checking. 
The matrix element of the uj — > Tr+Tr^Tr" decay is propotional to: 

{Ml-9mir ^' 

where: 

Pt^+ - momentum vector of 7r+ in the tt+tt^tt^ centre of mass, 
p^- - momentum vector of tt^ in the tt+tt^tt^ centre of mass, 
AfsTT - tt+tt^tt'^ invariant mass, 
- TT meson mass. 

If the distribution of A for background events is sufficiently uniform, then 
by weigthing events having A < 0.5 with -1 and those having A > 0.5 with +1 
we effectively subtract background. But the efficient number of signal events 
filtered this way is twice smaller than that in the initial sample. 

Both methods resulted in statistically compatible distributions of variables 
characterizing the wtt^tt" decay and production. The results of the subtraction 
procedure by the former method is presented in fig. ^. Two slopes are clearly 
seen in the momentum transfer squared distribution (t' = t — tmin) in fig- |^b. It 
is useful to divide the experimental event sample into two parts in the variable t': 
type LT ( low t' ) -t' < 0.08 GeF^ and HT ( high t' ) 0.08 < -t' < 0.7 GeV^. 
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Figure 2: a) Mass of the 7r+27r~27r°-systein with -t' < 0.7 GeV"^ from the 
band 0.7577 < M^+^-^o < 0.8077 GeF with superimposed distribution from 
the control bands, b) Momentum transfer t' for reaction n~Be — > u)Tr~Tr^Be . 
Mass of the Lvn~TT° system for: LT - c) and HT - d). 
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The p-meson peak dominates in the 7r~7r° invariant mass spectrum (fig. 
The 61(1235) isobar of the WTr-system is more evident in the high mass region 
of TT~^2-K~2-K^ (fig. ||b). The LO-K~n° invariant mass spectrum of the LT type 
(fig. ||c) has its maximum at 1.72 GeV and a weakly pronounced shoulder in the 
region around 2 GeV. The distribution of HT type events (fig. |^) is wider and 
has a shoulder in the region of 2 GeV. A peak cooresponding to the 02(1320) 
meson production is clearly seen in the low mass region. The PWA method 
was applied to study this complicated structure. Further the event sample, 
where the mass of at least one of the tt+tt^tt^ combination satisfies the relation 
\M^+^-^o - 0.78271 < 0.060 GeV", were analysed. 
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Figure 3: 
2.15 GeV 
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a) Mass of tt-tt*' with -t' < O.OSGeV"^, 1.65 GeV < M{-k+2t:-2t:^) < 
b) Mass of 7r+27r-7rO with -t' < O.OSGeV^, 2.15 GeV < 
-27r") < 2.45 Gey. 



4 Procedure of the partial wave analysis 

The goal of this method is the determination of intensities and relative phases 
of different states as a function of M^^-^o and t'. This method is based on the 
approach realised in Illinois PWA program |ll| . 

The full amplitude of each decay channel is written in the form of the product 
of production and decay amplitudes. The produced system decay process is 
assumed to be a chain of sequential two body decays into the isobar and a 
bachelor meson. 

The wTT^TT^-system state is described at fixed M(^^-^o and t' by the following 
set of quantum numbers: the total angular momentum J of unr^Tr^ system, the 
parity P, the projection of total angular momentum onto the Gottffried-Jackson 
reference frame z-axis M, the isobar orbital momentum in the LUTr~Tr^ decay L 
and the total spin S. It is convinient to change the basis of states having 
M = ±1 to the basis of their linear combinations, which are the eigenvalues 
of a reflection operation in the production plane. The quantum number which 
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characterises the behaviour of the system with respect to reflecton is called a 
naturality, 77 = ±1. As a result any partial wave is described by the following 
set of numbers: M'^ L[S]{isobar) with M > 0. 

A decay amplitude into the state with a definite orbital angular momen- 
tum L and a total spin S was written in the frame of a non-relativistic tensor 
Zemach formalism The isobar mass dependence of the decay amplitude 

was assumed in the relativistic Breit-Wigner form with the Blatt-Weisskopf 
form-factors 

Integrals of the amplitude products were calculated by the Monte-Carlo 
method using the set-up model. 

The selected events were divided into two parts in t': the LT and HT, each 
of which, in turn, was divided into the 50AIeV and 100 MeV bins in A^s^, 
beginning from — 1.2 GeV. The process of the beam diffraction dominates 
in the region of small momentum transfer, and this fact dictates the selected 
wave set with Af = 0+ and = 0", 1+, 2" The waves with AF' = 1+ and 

= 2"*", 4+ become dominant in the region of high momentum transfer. This 
is proved by the 02(1320) peak in the Ms^ spectrum for events from the HT 
sample. These facts determined the wave set used in the PWA. The following 
isobars were included: the in the 7r~7r° chanel, the bi, and pi, and /93(1690) 
in the ojtt chanel. The 1~1+S'(65) wave with indefinite isospin was included 
for the description of the visible 6° prevailing over 6^. The incoherent wave 
FLAT, uniformly distributed over the Stt phase space with one of the 7r+7r~7r'' 
combination from the oj region, was included to take the non LLnr~Tr'^ background 
into account. 

The resultant wave set looked as following: 
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The real and imaginary parts of the density matrix arc the parameters to be 
defined by the extended maximum likelihood method . The fits with density 
matrix of different ranks showed that the matrix of the rank 3 is sufficient for the 
experimental data description. The quality of description is illustrated in fig ^, 
where some angular and invariant mass distributions for events from one of the 
M^^-^o interval are plotted. The presented wave set was found to describe the 
data sufficiently well up to the mass M^^j^-^^o < 2.1 GeV. The background model 
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used in the PWA turned out to be good and the wave FLAT was smooth and 
its intensity was equal to the background in the uj peak region. The appUcation 
of another PWA method with the background subtraction by the variable A, 
mentioned above, proved the results of the PWA without subtraction. The 
results of the PWA of the bQMeV and \QQMeV binned data were consistent. 




cos(iJ) cos(iy) cos(i5) 

9) h) I) 

Figure 4: Comparison of the experimental data with events generated by 
the Monte-Carlo with the PWA model for -t' < 0.08 GeV^, 1.65 GeV < 
M(7r+27r-27r°) < 2.15 Ge^. a) and g) 7r+7r-27r", b) and h)7r+27r-7rO, c) and i) 
7r~7r''. d) Mass of 7r+7r~7r". f) Azimuthal angle of the 7r~-meson in the 7r~7r° 
centre of mass helicity reference frame. 



5 The results of the patial wave analysis 

The results of the PWA are presented in fig. || and ^. 
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Figure 5: The PWA results for events with —t' < O.OSGeV^. Wave intensities 
of: a) 0-0+Pl(/9), b) 2-0+P2{p), c) combined 1+0+ . Combined intensity of 
the l+0+P(6i) and l+0+Sl{p) is drawn in squares, d) Phase of the 0-0+Pl(p) 
wave with respect to the l+0+P(6i) wave. 
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J^M'' = 0~0+ : A peak in the region of 1.74^61^ with the flat background 
dominates in the wave 0-0+P(p) in the LT data sample (fig. |a). The same 
structure is less visible in the 0~0+P(p) wave intensity for events from the HT 
sample. The phase of these waves with respect to the smooth waves l"'"0"'"P(6i), 
l~^0'^S{pi) exhibits resonant behaviour ( fig. ||d). This fact makes it possible 
to assume the existence of a resonant state in the 0~0"'" wave. 

The resonance parameters were determined by the fit of the wave intensity 
with the incoherent sum of a relativistic Breit-Wigner function and a cubic 
polynomial. The parameter estimates are robust to the fit with a coherent 
sum and to the changes of the fit region. The mass was found to be M = 
1.737 ± 0.005 ± 0.015 GeV and the width T = 0.259 ± 0.019 ± 0.06 GeV. The 
number of resonance decays observed in the current experiment was estimated 
as Af(7r-(1740) cop') = (230 ± 12 ± 82) • 10^. The mass and width of this 
resonance differ from that for 7r(1800). The reason for this is unknown. It is 
possible to suggest the existence of two objects of difi'erent nature: a hybrid 
7r(1800) and a 3^5*0 qq state decaying into cup. 

= 1+0"*" : The intensity of the sum of all 1+0+ waves has its maximum 
at 1.75 Gey. Some structure can be seen at ^ 1.5 GeV^ (fig. ||c). This total 
wave is a result of the combination of strongly interfering waves: l+0+P(6i), 
l+0+S'l(p), l+0+5(/9i), l+O+Pi^p) and l+0+i:'2(p). The data do not allow 
one to state that the phase of some wave has a resonant behaviour in the region 
of - 1.5 GeF or 1.75 GeF. 

J^M^ — 2~0+ : A clear peak is observed in the region of mass 1.67 GeV with 
width ~ 0.2 GeV" (fig. |b). The resonant phase behaviour of the 2-0+Pl{p) 
and 2~0+P2(/9) waves with respect to the l+0+P(&i) wave makes it possible 
to conclude the observation of a resonance decaying into uip. This signal is 
nearly absent in the bin system. The resonance parameters of the 2~0+P2(p) 
peak were estimated in the same way as for 7r(1740): M = 1.687 ± 0.009 ± 
0.015 GeT^ and T = 0.168 ± 0.043 ± 0.053 GeT/. We identify this phenomenon as 
the decay of 7r2(1670) into ujp. 

The partial branching ratio was found by normalization to the decay 7r2(1670) 
/2(1270)7r 1^, observed in the current experiment decay: 

Br(7r2(1670)" ^ Lup-S = 2) = 0.019 ± 0.004 ± 0.010 
Pr(7r2(1670)- -> cup-S = 1) = 0.009 ± 0.002 ± 0.003 

The error of the normalization to the /27r decay is sufficiently large due to the 
Deck effect [|5|. The upper limits of the partial branching ratios of the decays 
7r2(1670) — > Trpi.pi — s- ujt: and 7r2(1670) — s- bin at the 2a confidence level were 
set: 

Pr(7r2(1670) pin) < 0.0036, 
Br(7r2(1670) ^ bin) < 0.0019. 
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It is worth mentioning that a relatively large probability of the 7r2(1670) — > 
ujp decay is in agreement with the model dependent calculations of 

J^ikf = 3+0+ : A peak at - 2 GeV of the width - 0.35 GeV is clearly 
seen in the 3+0+S'(/93(1690)) wave for events with low f (LT) (fig. |a). The 
resonant phase motion was not observed. Such wave behaviour can be attributed 
to the Deck effect process, characterized by the near threshold bump in the 5'- 
wave. A 3+0+ state of the up system is 3 times smaller at the intensity 
maximum. 

J^M^ — 1+1+ : This wave in the ujp^ is non-significant and does not have 
narrow structures of the width T 0.2 GeV. 

— 2+1+ : This wave dominates at high t'iWT sample). The intensive 
production and decay of the 02(1320) is the main process at low masses (fig. ^). 
Due to the strong interefence of the &i7r and ujp~ states in the region of their 
thresholds only the partial branching fraction of the a2(1320) W7r~7r° was 
determined relative to the well-known 02(1320) decay modes into rjir and p7r. 
This probability was found to be Br(a2(1320) wtt-tt") = (5 ± 1)%. There 
exists some arbitrariness in the question of what should be called the 02(1320) 
due to the complicated form of the ujtt~tt^ spectrum in the 2+1+ state. For 
possible misleading we define the 02(1320) partial width as that of a Breit- 
Wigner function with the S-wave oj'k~'k^ background. This decay was previously 
observed in the low statistics experiments . 

The nature of the 1.7 GeV mass structure is unknown. The simple suggestion 
is that this peak is due to the existence of a resonance with a mass M w 1.6 GeV 
and a width F « 0.3 GeF. The radial excitation of the 02(1320) meson should 
exist in this region. Nevertheless, it is worth saying that even in the absence of 
such a resonance the mass specrtum of the 2+1+5' tJ7r~7r° wave should have a 



complicated two bump shape due to the opening of ujp chanel |17 . Sufficiently 



good description of this two bump spectrum was obtained in the model of one 
resonance. These possibilities are different in the phase motion. The phase 
motion was measured in this experiment with large uncertainty and did not 
allow us to make a definite conclusion. The fraction of the 61 tt production in 
the same state in the band 1.6 < Ms^ < 2 GeV is negligible. 

jP M"^ — 4+I+ : There is a peak in the mass spectrum of top'' system with 
the spin 2 and orbital momentum 2 at w 2 GeV (fig. ^c). The phase shows 
resonant behaviour. This signal can be identified as 04 (2040). The 04 (2040) 
parameters were estimated by the mass spectrum fit with the incoherent sum 
of a Z)-wave relativistic Breit-Wigncr functon and a polynomial background. 
They were found to be: M = 1.944 ±0.008 ±0.050 Gel/ and F = 0.324 ± 0.026 ± 
0.075 GeF. The t' dependence of O4(2040) is identical to that of O2(1320). 

J^Af? ^ 1-1+ : This exotic wave has quantum numbers which are forbidden 
for qq states. The intensity of the l"l+S'(6i) wave is very small for the LT 
sample events. That for high t' shows a wide bump with the maximum at 
M - 1.6 ^ 1.7 GeV. The signal in the intensity of the J^A/'' = 1"!+ wave of 
the Ty'TT-system with close parameters was observed earlier in |l^ . All these facts 
indicate to the possibility of the existence of a wide resonance (F ^ 0.5 GeV) 
with mass M « 1.6 GeV . The maximum intensity of this wave does not exceed 
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15% of the 2"'"1+S'2(p) wave intensity (fig. ^). Similar waves in the up are 
also small. 

6 Conclusions 

The PWA of the W7r~7r°-system produced in the reaction ir^Be ii^2ti^2ti'^ Be 
with beam momentum = 37 GeV was performed. 

A resonance structure with M = 1.737 ± 0.005 ± 0.015 GeV and T = 0.259 ± 
0.019 ± 0.06 GeV was observed in the 0-0+ wave. These parameters 

differ from that of 7r(1800) seen in the 27r"7r+ system. This fact may point out 
to the existence of two different objects with close masses. 

The 02(1320)- wtt-tt" decay in the wave J^Af = 2+1+ was found to 
have relative probability i3r(a2(1320)^ ujtt^it^) = (5 ± 1)%. 

A wide bump of unknown nature at M « 1.7 GeV was observed in this wave. 

The 7r2(1670) — > Lop decay was observed in the J^iW = 2-0+ wave. The 
partial decay fraction was found to be: 

Br(7r2(1670)" ^ up-) = 0.027 ± 0.004 ± 0.01 

The 04(2040) up decay was observed in J^AI^ = 4+1+ wave. We mea- 
sured the 04(2040) resonance parameters : M = 1.944±0.008±0.050 GeV and V = 
0.324 ±0.026 ± 0.075 GeV. 

This research was supported in part by Russian Foundation for Basic Re- 
searches under grant no.98-02-16392 and INTAS-RFBR grant no.95-0267. 
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